GRADED ORBITAL OCCUPATION NEAR INTERFACES IN A La 2 Ni0 4 - La 2 Cu0 4 
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X-ray absorption spectroscopy and resonant soft x-ray reflectivity show a non-uniform distribution 
of oxygen holes in a La2Ni04 — La2Cu04 (LNO-LCO) superlattice, with excess holes concentrated 
in the LNO layers. Weak ferromagnetism with T c = 160 K suggests a coordinated tilting of NiOe 
octahedra, similar to that of bulk LNO. Ni d 3z 2_ r 2 orbitals within the LNO layers have a spatially 
variable occupation. This variation of the Ni valence near LNO-LCO interfaces is observed with 
resonant soft x-ray reflectivity at the Ni and Cu L edges, at a reflection suppressed by the symmetry 
of the structure, and is possible through graded doping with holes, due to oxygen interstitials taken 
up preferentially by inner LNO layers. Since the density of oxygen atoms in the structure can be 
smoothly varied with standard procedures, this orbital occupation, robust up to at least 280 K, is 
tunable. 
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I. INTRODUCTION 

Orbital states of transition metal (TM) oxides present 
the opportunity of adjusting material properties to a spe- 
cific purpose. In practice, spin and orbital reconstruc- 
tions by charge transfer have been observed near the 
LCMO-YBCO interfaces^, which locally have the struc- 
ture of 113 oxides between Cu and Mn atomic planes. 
Other layered structures have been proposed, where 
strain, electronic confinement or correlations remove the 
bulk orbital degeneracy^, leaving only one orbital avail- 
able for transport, similar to the Cu d x i_ y i orbital 
in high-temperature superconductors. The LaNi03 — 
LaA103 superlattices^— where the degeneracy of the e g 
orbitals in LaNi03 is removed, is an example of using 
these ideas for orbital reconstruction and control of cor- 
related phases of TM oxides with the 113 structure. 

It is intriguing to ask whether these ideas can be ex- 
tended to 214 materials. In some bulk 113 structures 
(e.g. SrTiOa) the TM ion moves away from the center 
of the oxygen octahedron and the orbital degeneracy is 
lifted by the crystal-field splitting. Strain can determine 
the relative order of thin film 113 e g orbitals if the or- 
bitals are degenerate in the bulkji For instance, when 
the apical oxygen atoms are moved farther along the c- 
axis, the d 3z 2_ r 2 orbital becomes lower in energy than 
the d x 2_ y 2 orbital. For TM ions compounds with the 
214 structure, the degeneracy of e g orbitals is similarly 
lifted by the crystal-field splitting. In this case, the order 
of e g orbitals is not modified by strain because the orbital 
energy change associated with it (~ 0.1 eV) is relatively 
small compared to the crystal- field splitting (~ 1 eV). In 
addition, in contrast to 113 structures, additional atomic 
layers are present between TM layers in 214 compounds. 
This raises the question of whether 214 heterostructure 
interfaces will affect the occupation of TM orbitals away 
from bulk properties sufficiently strongly, through con- 
finement effects, charge transfer or in other ways. 

Here, we investigate the effect of interfaces on oxygen 
dopant distribution and TM orbital occupations in a su- 



perlattice made of La2Ni04 (LNO) and La2Cu04 (LCO) 
layers. We observed that oxygen atoms are doping pref- 
erentially the LNO layers, changing the TM orbital oc- 
cupation. More importantly, they can allow a reflection 
suppressed by symmetry, with the dopant-induced scat- 
tering separated in momentum from other charge scat- 
tering. A spatially variable Ni valence and d-orbital oc- 
cupation within the LNO layers is observed in this way 
with resonant soft x-ray scattering. We therefore found 
that superlattice interfaces can significantly affect the oc- 
cupation of the TM orbitals in 214 oxides. 

This corroborates the observation that oxygen content 
can determine the transport properties of 113 TM oxide 
heterostructures. For instance, LaA103 films on SrTi03 
(STO) substrates are insulating or metallic, depending on 
the number of oxygen vacancies^ Our results show that, 
to understand the properties of oxide heterostructures, it 
might be necessary to consider the oxygen dopant distri- 
bution, not only the average content. Since the oxygen 
interstitial distribution is relatively easily modified with 
standard procedures, this orbital occupation in 214 LNO- 
LCO superlattices is tunable. 



II. EXPERIMENTS 
A. Superlattice structure 

The superlattice, made of LNO and LCO layers on a 
STO substrate, was grown by molecular beam epitaxy 
at Brookhaven National Laboratory. Atomic force mi- 
croscopy (AFM) measurements on a Dimension 3100 in- 
strument showed good superlattice surface quality with 
an average rms roughness over a 5 (im x 5 fim area of 
0.53 nm. The sample has La in both layers; therefore 
it is not possible to use scattering near the La edge to 
characterize the interface roughness. 9 

Off-specular hard x-ray scattering measurements (not 
shown) confirmed that the superlattice is fully strained 
in the a — b plane to the STO substrate with a s = b s = 
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FIG. 1: (color online) (a) STEM high angle annular dark field 
image of a LNO-LCO superlattice showing epitaxial growth. 
The right panel expands the image area enclosed by the rect- 
angle. The bright features in the image are La atoms in the 
LaO planes, represented by dashed lines between the TM — 06 
octahedra. The oxygen atoms are not visible. The experi- 
mental conditions were 1.3 A probe size at 200 kV, 28 mrad 
convergence angle, and 64 — 341 mrad collection angle, (b) 
Wide angle hard x-ray measurements made at Advanced Pho- 
ton Source at Argonne National Laboratory (7.5 keV) and 
with an X' Pert diffractometer (8.05 keV) showing superlat- 
tice peaks (SL) near the main reflections (M) indexed with 
respect to the measured superperiod. Substrate peaks (STO) 
are also visible. 



3.915±0.035A. A comparison to measurements on bulk 
LNO and LCO single crystals*!^ indicates that the LNO 
layers are under tensile stresses of 1.5 % and 0.6 % along 
the a and b directions of the STO substrate, while the 
LCO layers are under tensile stresses of 3.4 % and 2.7 % 
along the a and b directions. 

Hard x-ray reflectivity (HXR) scans (Fig. lb) , made at 
beamline ID6 at the Advanced Photon Source at Argonne 
National Laboratory and using a laboratory X'Pert x-ray 
source, showed pronounced superlattice peaks. The aver- 
age monolayer thickness was 6.470 ± 0.035 A, defined as 
1 "molecular layer" (ML), and corresponding to half the 
average of LNO and LCO unit cells c-axis parameters. 
The number of ML in a superperiod was determined by 
measuring the ratio of momentum transfers for superlat- 



tice reflections satellites (SL) to momentum transfers of 
Bragg peaks (M in Fig. lb) corresponding to 1 ML. The 
average of ID6 and X'Pert measurements showed a super- 
lattice period of (7.45 ± 0.08) ML and (7.42 ± 0.14) ML, 
respectively, indicating an approximately equal mixture 
of 7 ML and 8 ML superperiods. 

The average ML thickness, number of MLs in a su- 
perperiod and number of superperiods (12 deposition re- 
peats), give a total superlattice thickness of 577 ± 10 A, 
consistent with the result obtained from thickness oscil- 
lations (Fig. 2) of 570 ± 20 A. 

The superlattice deposition sequence had LNO and 
LCO layers of nominally equal thicknesses. Indeed, HXR 
measurements showed a step- like feature at L = 2, with 
the scattering momentum Q — 2irL/c indexed with re- 
spect to the superlattice superperiod c for low-Q range, 
on a large sloping background (ID6 curve in Fig. 2). 
This arises from the slightly different c-axis parameters 
of LNO and LCO unit cells. Scattering at the L = 2 
superlattice reflection is suppressed for equal LNO and 
LCO layer thicknesses. Thus, the thicknesses of the LNO 
and LCO layers will be considered equal in the following 
analysis. However, the L = 2 peak, suppressed for non- 
resonant HXR scattering, is observed for resonant x-ray 
scattering at the Ni edge, but not at the Cu edge. This, 
regardless of any irregularities of the structure, highlights 
a difference in the initial and intermediate state distri- 
butions of resonant scattering at the Ni edge, that is a 
difference between the distributions of atomic cores and 
Ni valence electronic states. In Section II. E, it will be 
shown how this can happen at the Ni edge (Fig. 2). 



B. Superlattice magnetism 

Undoped bulk LNO and LCO are antiferromagnetic 
insulators with two and one unoccupied TM orbitals and 
nominal 5=1 and 5=1/2 spins on the Ni 2+ and Cu 2+ 
ions, respectively. Bulk La2NiO4.03 shows a weak ferro- 
magnetism below 65 K with a net magnetic moment ori- 
ented along the c-axis of approximately 0.04 ^b/Ni ion 
at 4.2 K. 1 - The ferromagnetism is not the consequence 
of oxygen doping; the NiOg octahedra in undoped bulk 
LNO tilt in the low-temperature tetragonal phase, mis- 
aligning the spins and making the material weakly ferro- 
magnetic. 13 In doped LNO the ferromagnetic phase ap- 
pears below - 150 - 200 Ki 3 . 

SQUID measurements of the superlattice, done with 
a Quantum Design instrument, showed a ferromag- 
netic moment below T = 160 K (Fig. 3), with < 
0.05 ^e/(Ni ion), consistent with the magnitude of bulk 
LNO ferromagnetism. The SQUID measurements also 
did not show a superconducting anomaly in the mag- 
netic moment down to 10 K, which sets an upper limit 
on hole doping of the LCO layers of 0.06?^ 

The consistency of the magnitude of the moment and 
temperature of the transition with that of bulk LNO sug- 
gest that the superlattice ferromagnetism is likely con- 
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FIG. 2: (color online) L scans of the low momentum range 
for hard and soft x-ray scattering. HXR measurements show 
strong thickness oscillations, with a peak at L — 1 and a step 
at L = 2 (black line) for ID6 data. The peaks at L=3.8, 
4.6, and 5.7 are stray reflections. Soft x-ray scattering shows 
pronounced superlattice (SL) peaks for L — 1 at O K, Ni L 
and Cu L edges. The Ni L = 1 peak is shifted because of 
refraction. The L=2 reflection is much larger at the Ni edge 
compared to the Cu edge. This observation is confirmed by 
detailed energy-resolved measurements (Fig. 6) which are dis- 
cussed below. In contrast, the absence of the L—2 reflection 
in measurements at the oxygen edge is caused by the energy 
used (see Fig. 5 for a complete measurement). Smaller dif- 
ferences can be observed between the measurements at the 
Cu and Ni edges, for instance a slight indication of a peak at 
L = 4 at the Cu edge. The small difference in intensity at 
L = 4 can be explained using the different polarizations of Ni 
and Cu d-orbitals, unlike the opposite and larger difference 
at L — 2 (Sec. II E). The intensity of SL reflections at L — 3 
and L = 4 is small, consistent with the calculated dependence 
of SL structure factors amplitude on L. The L — 5 and L — 6 
SL peaks are larger because of the proximity to the Bragg 
peak corresponding to I ML (off a;-axis in this figure). 



fined to the LNO layers, where the magnetic moment 
has a similar origin to that of bulk LNO. The necessary 
corresponding modifications to the magnitude of Ni spins 
near the interfaces from different Ni orbital occupations, 
discussed later and shown schematically in Fig. 8b, do 
not appear to result in large additional ferromagnetic mo- 
ments. In this case, this is likely because exchange cor- 
relations decay spatially very fast and the net magnetic 
moment of Ni02 and Cu02 layers in bulk compounds is 
small or zero, respectively. The temperature variation of 
the SL ferromagnetism will invalidate a candidate expla- 
nation for the difference at L=2 in the resonant soft x-ray 
scattering at the Ni and the Cu edges (Sec. II E). We 
will use the smallness of the superlattice ferromagnetic 



moment in Section III. A. 



C. Hole distribution 

Oxygen doping bulk LNO and LCO places oxygen 
interstitials between the LaO layers of the structured 
The electronic structures and orbital orientations of bulk 
LNO and LCO, studied with x-ray absorption spec- 
troscopy (XAS) experiments on bulk nickelate a 13 ' 16 ' 17 
and cupratesi*, showed doping-dependent upper Hub- 
bard bands (UHB) and mobile carrier peaks (MCP). 

XAS measurements showed a non-uniform distribution 
of holes in the SL and the favorable conditions for a large 
Ni valence modulation. We used XAS in total fluores- 
cence yield (TFY) and total electron yield (TEY) modes 
to characterize the hole distribution in the LNO-LCO su- 
perlattice (Fig. 4a). Measurements were made at beam- 
line X1B at the National Synchrotron Light Source at 
Brookhaven National Laboratory. To determine which 
peak belongs to which layer, the short probing depth 
of TEY compared to TFY was used to isolate the sig- 
nal from the shallower layer of the superperiod, which is 
known from the deposition sequence to be LNO. For en- 
ergies in the range of a few hundred eV, as in our case, the 
TEY probing depth is less than 100 A. In contrast, the 
penetration depth of TFY at these energies is thousands 
of A, effectively probing the entire superlattice. There- 
fore, the peaks that are present in TFY but not TEY 
originate from the deeper layer of the superperiod pair, 
which is the LCO layer. 

Fig. 4a shows TFY and TEY measurements on the 
superlattice at room temperature. The lowest energy 
peak that originates in the LCO layers aligns with the 
upper-Hubbard band (UHB) peak of bulk LCO^. while 
the lowest-energy peak in the LNO layers aligns with 
the mobile carrier peak (MCP) of bulk LNOJ&i 7 This 
suggests that the LCO layers of the superlattice are un- 
doped, consistent with the magnetization measurements 
in Section II. B, while the LNO layers are doped by excess 
oxygen. The LCO MCP peak begins to be clearly present 
at x ~ 0.03 in bulk La2- x Sr x CuO,fi£; since no such peak 
is visible for the LNO-LCO superlattice, it must be that 
xlco < 0.03. Similarly, the LNO UHB is present up to 
x ~ 0.2 in bulk La2_ x Sr x NiO,fi&; since no peak is seen 
in TFY or scattering other than LCO UHB in this en- 
ergy range, xlno > 0.2 or oxygen doping 5 > 10 %. We 
conclude that the doped holes are localized in the LNO 
layers of the superlattice. 

Even with a large number of holes in the LNO layers, 
the superlattice remains insulating (Fig. 3). Assuming 
that holes are shared equally between Ni and O atoms, 
xlno > 0.2 implies that the average Ni valence increases 
from 2+ in undoped LNO to larger than 2.1+ in oxygen- 
doped LNO. It can, in principle, accommodate a substan- 
tial modulation of the Ni valence within the LNO layers. 
This modulation is observed by scattering at the Ni edge 
(Section II.E). 
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FIG. 3: (color online) Magnetic moment measurements show 
a ferromagnetic transition at 160 K. A diamagnetic contribu- 
tion has been subtracted from the data. The applied magnetic 
field was 100 Oe. Resistivity measurements show an insulat- 
ing behavior down to 38 K with a small anomaly in slope at 
58 K. 



D. Scattering at the oxygen edge 

Resonant soft x-ray scattering (RSXS) can probe bulk 
charge order.— ~— Its energy and momentum selectiv- 
ity makes it a useful technique for also probing the 
spatial distributions of orbitals of different energies in 
nanometer-sized structures. Measurements on the LNO- 
LCO superlattice at the oxygen K, Ni and Cu L edges (see 
Fig. 2 for selected resonant energies) were made at beam- 
line X1B at the National Synchrotron Light Source in an 
ultra-high vacuum diffractometer with a 7r-polarized inci- 
dent beam. Scattering measurements were made at 90 K, 
unless specified otherwise. The beam was focused to an 
approximately 1 mm x 1 mm spot, which is much smaller 
than the sample size. 

RSXS at the oxygen K edge probes the spatial distribu- 
tion of the occupation of oxygen p-states and illustrates 
one reason for the L dependence of resonant intensities 
of SL reflections. Measurements were made in specu- 
lar geometry with a variable energy E and momentum 
perpendicular to the surface Q = 2ttL/c, with reflec- 
tions indexed with respect to the superperiod c for low-Q 
range (for other examples of RSXS on superlattices, see 
Figs. 4b and 5 show reflectivity measurements at 
the oxygen edge for L=l and L=2, respectively. Refrac- 
tion effects and small superlattice imperfections slightly 
shift the peaks from integer values. Peaks at 529 eV 
and 531.5 eV are visible for the momentum held fixed 
at L=l (Fig. 4b), indicating a spatial modulation of 
the same LNO MCP and LCO UHB states that are visi- 
ble at 528.5 eV and 530.5 eV in FY measurements (Fig. 
4a). A third peak is present for L=2 at 533.5 eV (Fig. 
5); this peak will be addressed in Section III.B. The ra- 
tio of the scattering intensities for LNO MCP and LCO 
UHB energies is different for L=l and L=2. Specifically, 
Ilno,mcp/Ilco,uhb = 1.1 for L=l and 0.4 for L=2. 

The intensity of a superlattice reflection in the single- 
scattering approximation is given by I — A\e*fSL^i\ 2 , 



where 

s£'=£<*r>/;w 2,rWc (i) 

/ , n 

is the structure factor, (t™) is the average over a layer of 
the occupation factor i™ (equal to or 1), which identi- 
fies the valence and element n to use in Eq. 1 for each 
site in layer zi (e.g. Ni 2+ or Ni 3+ ), and f^(ui) is the 
atomic form factor for valence and element n, with the 
indices i and j accounting for X-ray birefringence. Defin- 
ing the variables {tf 1 } in Eq. 1 separates the spatial (I) 
and energy (to) dependence of the scattering factor. This 
is justified by the very weak dependence of atomic form 
factors on momentum for soft x-ray energies. In compar- 
ing the intensities of LNO MCP and LCO UHB features, 
we neglect self-absorption effects and the change with ge- 
ometry of A (e.g., footprint effects) since the scattering 
geometries and energies are very close. The variation 
with L of the intensity ratio Ilno,mcp/Ilco,uhb can, in 
principle, be due to either the polarization dependence 
of the form factors f% or the spatial distribution of the 
occupation factors t™. 

Bulk LCO has hole orbitals oriented in the a — b 
planed, while the holes in bulk LNO have a substantial 
component oriented along the c-axi o 17 ' 24 ' 27 . The varia- 
tion of the intensity ratio Ilno,mcp/Ilco,uhb with L can 
be qualitatively explained by assuming a polarization of 
hole orbitals in the superlattice LCO and LNO layers 
that resembles that of bulk LCO and LNO. Specifically, 
the scattering intensity of a horizontally-polarized beam 
off intermediate states oriented along a — b and c-axis is 
proportional to sin 4 (#) and cos 4 (8) respectively, where 
9 is the angle between the incident beam direction and 
the sample surface. The intensity for the LCO state will 
be proportional to sin 4 (#). In contrast, for the LNO in- 
termediate state, with a substantial component oriented 
along c-axis, the scattering intensity has a term propor- 
tional to cos 4 ((9) and can be expected to increase slower 
or even decrease as 8, and therefore L, increases. Then, 
to explain the variation with L of Ilno,mcp/Ilco,uhb 
observed in the measurements, it is sufficient to assume, 
based on the relatively small influence of strain in 214 
structures (Section I), that the hole p-orbital orientations 
in LNO and LCO layers are similar to those in bulk crys- 
tals. 

The oxygen p orbitals hybridize with the TM orbitals 
inside the TM — Oq octahedra (Fig. la). Therefore, the 
polarization of the oxygen form factors / y is related to 
that of TM form factors and a similar variation of the 
scattering intensity with L at the TM edges is expected. 
However, in contrast to the case of oxygen edge LNO 
MCP and LCO UHB states, in the next section it will 
be shown that the evolution with L of the scattering in- 
tensities for Ni and Cu orbitals cannot be explained with 
polarization of form factors alone. 
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FIG. 4: (color online) (a) Electron and fluorescence yield and 
the correspondence of the peaks to the superlattice layers. 
The LCO UHB peak disappeared completely for an LNO- 
SCO superlattice grown on an STO substrate, showing that 
there is no STO substrate or LNO layer contribution at the 
LCO UHB energy. From detailed measurements (not shown) 
on LNO-SCO superlattices LNO MCP and LCO MCP have 
different energies, (b) RSXS intensity (logarithmic scale) at 
90 K at the O edge shows two peaks corresponding to the first 
two peaks in TFY, with the lower (529 eV) and higher (531.5 
eV) energy peaks from the LNO and LCO layers, respectively. 
The different L profiles are due to an interference effect and 
different layer c-axis parameters. 



E. Scattering at transition metal edges 

Resonant scattering at TM edges probes the d-states 
of Cu and Ni ions. L scans on resonance at these en- 
ergies (Fig. 2) show the expected peaks at L=l. The 
L=2 peak is not visible at the Cu L edge in Fig. 2; 



however, the L=2 peak is unexpectedly pronounced at 
the Ni edge. Detailed energy plots of this momentum 
region, shown in Fig. 6, confirm the large difference 
I(L = 2)/I(L = l)| Ni >> I(L = 2)/I(L = l)| Cu between 
the L=2 reflections at the two edges: the peak at the Ni 
edge is 100 % more intense than the background, com- 
pared to only 6 % at the Cu edge. This superlattice reflec- 
tion was observed for sample azimuthal orientations with 
incident beam polarization at 45° and along the tetrago- 
nal a — b axes (Fig. 7a). In contrast to the temperature 
dependence of the magnetic moment in Fig. 3, the in- 
tensity of the L—2 peak at the Ni edge does not change 
with temperature between 90 K and 280 K (Fig. 7b). 
Therefore, this peak is not related to ferromagnetism. 

Contrast the observation of the larger L=2 intensity 
at the Ni edge compared to the Cu edge with the oppo- 
site trend at the O edge for LNO MCP and LCO UHB 
states, in which the LNO state intensity is suppressed at 
larger L. Therefore, the difference in the intensity ratio 
I(L = 2)/I(L = 1)|ni,Cu f° r scattering at Ni and Cu edges 
cannot be explained using the difference in the orienta- 
tion of p-states in the LNO and LCO layers, noted at the 
oxygen edge and a similar polarization of d-state TM or- 
bitals hybridized with oxygen orbitals in their respective 
layers. The polarization of the form factor f%(co) cannot 
be the explanation; instead we must look to a difference 
in the occupation factors t™ for Ni and Cu d-orbitals. 

This difference originates in the different dependence 
of the Ni and Cu valence on doping. LNO resembles a 
Mott- Hubbard insulator more than LCO.— This is vis- 
ible in the small degree of contrast in TFY at the Cu 
edge with doping of bulk LCQ2& and the large contrast 
at the Ni edge with doping bulk LNO—. Therefore, the 
Ni ci-orbitals participate in the hole doping process more 
than the Cu d-orbitals. The Cu d-states occupation fac- 
tor tf follows the distribution of the Cu ions because their 
valence remains close to Cu 2+ . In contrast, the substan- 
tial oxygen doping in the LNO layers (Section II. C) and 
the strong coupling between hole doping and Ni d-states, 
empties additional Ni orbitals, thus raising the Ni ion 
valence toward Ni 3+ . 



Therefore, unlike the Cu case where oxygen doping is 
absent, the spatial distributions of both Ni atoms and 
oxygen interstitials determine the occupation of Ni in- 
states that is, the magnitude of tf. If the oxygen dopants 
are not uniformly distributed within the LNO layers, a 
variable Ni valence results, i.e. spatial distributions tf 
for the doped orbitals different from those of the Ni ions. 
In this case, as pointed out in Section II. A, since the 
intermediate d-states do not follow the distribution of the 
underlying structure represented by the core atoms, the 
scattering at L=2 is not suppressed. In the next section, 
we will use this observation and Eq. 1 to determine the 
orbital occupancies i™ in the LNO layers. 
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FIG. 5: (color online) Intensity of the L—2 reflection at the 
O edge at 90 K (logarithmic scale). The three peaks are 
(in order of increasing energy): LNO holes at E=529.5 eV, 
LCO UHB holes at E=531.5 eV and a feature at E=533.5 
eV, discussed briefly in Section III.B. Unlike the scattering at 
L=l, the LNO MCP peak has substantially smaller intensity 
than the LCO UHB peak. 



III. DISCUSSION 

Our main finding does not depend on structural ir- 
regularities. The distributions of Cu and Ni atoms, and 
therefore of core p-states, are complementary (the layer 
occupations are {tf u ) + {t^ 1 ) = 1, for each layer I, where 
(t) is the occupation factor from Eq. 1), as supported by 
the observation of a crystalline SL structure with Cu and 
Ni atoms occupying their normal sites within the 214 unit 
cell. Complementary distributions have the same abso- 
lute value of the Fourier transform amplitude at L=2. 
Therefore, if one does not consider differences in the Cu 
and Ni d-state occupations, since the core p-states (the 
initial states of the RSXS transition) distributions are 
complementary, for any value of superlattice superperiod, 
specific model of the distribution of interface atoms (not 
necessarily Gaussian as we assume in Sec. Ill B) or size of 
interface roughness, the L=2 peak cannot be suppressed 
at the Cu edge without being also suppressed at the Ni 
edge. To explain the difference in scattering at L=2, it 
is necessary to consider differences in the (i-states (the 
intermediate states of the RSXS transition) that is, in 
the electronic properties of the LCO and LNO layers. 

The difference in the L=2 scattering at the Cu and Ni 
edges cannot be attributed to a different resonant form 
factor f l i(w) for the Cu and Ni edges. The energy u> and 
momentum L variables are separable in the soft X-ray 
scattering cross-section (Eq. 1); therefore, since Cu and 
Ni atoms scatter at L=l, considering only the different 
form factors /^(w) will imply that both scatter at L=2, 
contrary to observations. 

The conclusion that the Ni valence is not constant in 
the LNO layers is obtained from comparing the presence 



of L=2 scattering at the Ni edge with its absence at the 
Cu edge in the scattering model of Sec. Ill A. The mea- 
sured ratio I(L=2)/I(L=1) at the Ni edge and use of spe- 
cific SL structure models are necessary for a quantitative 
estimate (Sec. Ill B). 



A. Scattering model 

Calculating the energy dependence of Ni edge scatter- 
ing requires a knowledge of the energy and doping de- 
pendence of the scattering factors / ; IJ,Nl (w) in the LNO 
crystalline lattice. Band structure calculations of the Ni 
density of states variation with doping with LSDA+U 
models^ cannot reproduce the XAS experiments in the 
low-doping ranged. A detailed cluster analysis of en- 
ergy levels of Ni 2+ and Ni 3+ ions in an octahedral crys- 
tal field shows that the degeneracy of the x-ray transi- 
tions is lifted, which results in a very complex transition 
spectrum32 This is visible in Fig. 8a (top) for Ni 2+ and 
Ni 3+ ions in a cubic field. The scattering at L=2 at the 
Ni L3 and L2 edges is nevertheless made of only one peak 
(Figs. 7c inset and 8a). This suggests that a simplified 
scattering model is warranted. 

The scattering at the Ni edge can be analyzed using 
Eq. 1 and a simplified model. We define the functions 
fn(w) first. The XAS and L — 2 scattering energy pro- 
file (Fig. 8a) do not resemble the density of states of 
the intermediate unoccupied Ni d-states; for instance, 
each of the two peaks in the XAS at the Ni L3 edge, 
following immediately after the La edge at 853.5 eV in 
Fig. 8a, has both x 2 — y 2 and 3z 2 — r 2 character^ 
However, the state into which the electron is excited in 
the transition is well-defined and will be used as a la- 
bel. Therefore, the transitions at the Ni L3 edge are 
labeled by the intermediate d-states involved in the 2- 
step RSXS process and organized into four main groups 
based on the spin of this intermediate state (up or down) 
and its orbital symmetry (d 3z 2_ r 2 and d x 2_ y 2). For un- 
doped LNO (or a Ni 2+ valence), two of these states are 
unoccupied (d 3z 2_ r 2 ^ and d x 2_ y 2+). Doping toward a 
Ni 3+ valence is a complex process that cannot be repre- 
sented as sequential emptying of additional orbital levels. 
However, because of the isotropic nature of the holes in 
bulk LN O 17 ' 24 , it is expected that both the other two 
d-states (d 32 2_ r 2 ± and d x 2_ y 2^) will be emptied. In- 
deed, polarization-dependent measurements of bulk LNO 
at the Ni edge show that these states are approximately 
equally emptied for < x < 0.1 322, becoming available 
for the scattering process. As at the oxygen edge, it is as- 
sumed that the bulk / lJ can be used for the layers of the 
LNO-LCO superlattice since, as pointed out in Section I, 
strain will not change the relative order of d-orbital en- 
ergies in compounds with the 214 structure. Therefore, 
the atomic form factors for Ni ions in the layers are: 
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FIG. 6: (color online) (a) Scattering intensity (logarithmic 
scale) at La M 4 (E=853.25 eV) and Ni L 3 (E=856 eV) edges 
at L—2. The increase in the intensity over background at the 
Ni edge is 100 %. The features below 853 eV are thickness os- 
cillations. No off-resonant intensity tail is observed near L — 2 
at the Ni edge, showing that this reflection originates in an 
electronic state and not from, for instance, selective removal 
of Ni ions from the LNO layers, (b) Scattering intensity (log- 
arithmic scale) over a wide range at the Cu edge near L — 2 
for comparison. The energy range (14 eV) is chosen the same 
as in (a), (c) The area of interest at the Cu edge magnified, 
compared to TFY and TEY. A very small peak is visible for 
L=1.96, with a relative change in intensity over background 
of only 6 %. Because of the large sloping background an L 
scan is less useful in showing this small peak. 



4 + fl 



(2+)-!-(3+) 



(3) 



where f^+i^) 



is the form factor for Ni 2+ ions, 



t%3 
■'(24 



>( - 3+ ^(cj) is the term added with increasing the Ni 

valence to Ni 3+ , and the off-diagonal terms have been ne- 
glected because the ferromagnetic moment is small (Sec- 
tion II. B). These equations reflect the observation that 
in undoped LNO the d 3z 2_ r 2,^ and d x 2_ y 2,\ levels are 
empty and that doping empties in approximately equal 
amounts the orbitals of the opposite spin. 

The other term in Eq. 1 is (£"). To model the 
momentum dependence of scattering through the or- 
bital occupancies i™, we introduce the variables p 3+ and 
pf + = si — p 3+ , representing the layer-resolved average 
number of Ni atoms with nominal Ni 3+ and Ni 2+ valences 
in 1 ML over a a s x b s area, where s; is the distribution 
of Ni ion cores. For instance, the probability that the 
Ni ions have the nominal Ni 3+ valence is represented by 
< p 3+ < s;. From Eq. 3, the d 3z 2_ r 2 7 \. and d x 2_ y 2,\. 
levels are unoccupied with a distribution following that of 



' ' Thus, the average over each layer of the occupation 



factors for Ni ions are (i ; 



si, (t. 



si, 



„3+ 



<*r r ' l )=Pt + and (t? ^)=pj 

Then, the average Ni scattering factor in layer I de- 
pends on the variable number of unoccupied orbitals as: 



(2+)-K3+) 



(4) 



where f2 + (uj) has the distribution of Ni cores {s;}, 



to 



which spatially-dependent doping adds an additional 
„3+ „„ A f ij /^y Therefore, 



term proportional to p 3+ and fn+y 



<3+V 



the spatial symmetry of f, J ' '(w) is determined not only 
by the distribution of Ni core states {si}, but also by 
the intermediate d-state through the dopant distri- 
bution. The symmetry requiring the L—2 scattering to 
be small, noted in Sections II. A and II. E, is lifted in this 
way, by making the valences, and therefore the scattering 
factors, of Ni ions in outer and inner LNO layers different 
(Fig. 8b). This allows the L=2 reflection at the Ni edge. 
In contrast, because there is no change in the valence 
of the Cu ions (Section II. E), the Cu scattering factor 
u (lu) depends on structure only, through {1 — s;}, 
the complementary distribution to {s;}, and the L = 2 
reflection remains suppressed. 



B. Ni valence modulation 

We use this model to arrive at an estimate of the mag- 
nitude of Ni valence modulation in the LNO layers. Self- 
absorption is neglected because the calculated attenua- 
tion lengths are 1400 A and 1600 A at the Ni and Cu 
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an increase in the intensity above the edge*i, not seen 
in measurements on the LNO-LCO SL at the Cu edge 
(Fig. 6c). The slope in FY in Fig. 7c (inset) is due to 
a peak in the undulator spectrum. The relatively small 
self- absorption and the RSXS scattering geometry away 
from grazing incidence (the detector angles were 19.8° 
and 37.6° for L=l and L=2 reflections at the Ni edge) 
support the use of a single-scattering model for the anal- 
ysis of the RSXS measurements. 

The contribution to scattering from, for instance, far- 
away resonances and the relatively small difference be- 
tween the c-axis parameters of LNO and LCO is negligi- 
ble because it would give a non-resonant energy intensity 
tail at L=2 near Ni edge, which is not observed (Fig. 6a). 
S2, the L=2 structure factor in Eq. 1, will then be pro- 
portional to A = I J2iP^ +el4 '* z '^ c \> the L=2 component 
of the Fourier transform of the set 

In addition, S2 will be proportional to /(2+)^(3+), giv- 
ing the energy variation of L = 2 scattering. However, 
the L=2 scattering is clearly made of only one peak (Figs. 
7c and 8a) and the relatively small scattering angles sup- 
press the scattering intensity from the doped d x 2_ y 2,\. 
states, which suggests that we can neglect this scattering 
channel. The contributions to scattering from terms pro- 
portional to f xx and f vv are suppressed at small scatter- 
ing angles and the contribution from a small er-polarized 
beam component is neglected. Then, the only terms left 
from Eqs. 2 and 3 are those proportional to /J z 2 2 ^ 

and /J3,2-,2,l- 

The discussion in Sec. Ill A is independent of the 
specific details of the SL structure, given by the {sl} 
distribution. However, to obtain a numerical estimate, it 
is necessary to refer to the SL structure from our HXR 
measurements. The intensities of L=2 and L=l scatter- 
ing for an 8 ML superlattice at the Ni edge given by Eq. 
1 are: 



I(L = 2)\ m = A L=2 (Af 



COS 



?L= 2 ) (5) 



FIG. 7: (color online) (a) L scans at the Ni L3 edge show the 
L — 2 peak for sample azimuthal orientations with incident 
beam polarization at 45° and along the tetragonal a — b sixes, 
(b) The L = 2 peak at the Ni L3 edge persists up to at least 
280 K. As will become clear in Section III. A, this is due to 
the larger energy necessary to move the oxygen dopants in the 
structure compared to the thermal energy needed to destroy 
a spin order, (c) The L = 2 peak at the Ni L2 edge. The 
inset shows the energy profile of this peak (circles) compared 
to TFY (squares) and TEY (blue line). 



edges, respectively, compared to the total SL thickness 
of 575 A, and the measured inverse momentum widths of 
the L=l SL peaks, equal to 525 ± 50 A and 620 ± 50 A 
at the Ni and Cu edges respectively, are consistent with 
the total SL thickness of 575 ± 20 A. In addition, self- 
absorption has a clear signature in FY in the form of 



I(L = 1) 



Ni 



A L=1 (2.61) 2 \fr 



2 cos 4 ( 



,=i) (6) 



where | J2Li si^ Zllc \ = I Ef=i ^ Zljc \ = 2.61 and the 
smaller contribution to S\ proportional to A has been ne- 
glected. Using /|* a _ r3 4 ~ fdl 2 2 , t (from, for instance, 
the small XMLD "for Ni 24 
print effects (Al = i ^4l=2) 
scattering intensity is: 



I 3i 2_ r 2,t 

ions^ 2 -) and neglecting foot- 
the ratio of L~2 to L=l 



I(L 



I(L = 1) 



(#) 2 cos 4 ( 



'L=2 



(2.61) 2 cos 4 (6» L=1 ) 



(7) 



The scattering at L=2 is well separated in energy at the 
Ni and La edges (Fig. 6). However, this is not the case 
for L=l.— In practice, to obtain the value of the scatter- 
ing intensity for L=f at the Ni edge, we cut through the 
2-dimensional profile at the energy where the L=2 scat- 
tering has the maximum intensity. This is justified by the 
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FIG. 8: (color online) (a) Calculated XAS (offset vertically for 
clarity) using the software "Missing"~ for the Ni 3+ and Ni 2+ 
valences are in good agreement with Ref.—. The parameters 
used were lQDq = 1 eV, Ds = Dt — 0, 0.8 Slater integrals 
scaling, 0.2 eV Lorentzian and 0.1 eV Gaussian broadenings. 
Instead of the full tetragonal D^h symmetry, which further 
splits the final d-states, a cubic crystal field, with d-states of 
octahedral Oh symmetry, has been used for simplicity. The 
lower panel shows the measured electron yield (navy line), 
fluorescence yield (red squares) at La M 4 and Ni L 3 edges, 
and the fit of TFY (red line) with three peaks at the La and 
Ni edges (green) on a smooth step function (blue). The scat- 
tering at L=2 (black squares) peaks at the first Ni peak only. 
The lowest energy peak at 851 eV is a fringe effect, seen more 
clearly in Fig. 6a. (b) Sketch of the layer-resolved orbital oc- 
cupation where, for clarity, a LNO layer with 4 ML has been 
illustrated and only the rf3 Z 2_ r 2 orbitals have been included. 
Dashed lines represent the double LaO layers between which 
oxygen interstitials are located. Their number has been exag- 
gerated in this sketch and the tri-lobe drawing (not to scale) 
illustrates their hybridization with La ions in the two adja- 
cent LaO layers along axes at 90° with each other.— A small 
spin gradient is also shown at the interface. 



simple shape of the SL resonance profiles at the Ni edge. 
The small contribution from the La edge tail is neglected. 

From the measurements tttEtI = 0.0185, 6l=i = 9.9° 
and L=2 = 18.8°. Then, from Eq. 7, A = 0.77, or a 
change in the average Ni ion valence between the inner 
and outer layers of (p 3 + - p 3 ^ t ) = A/%/2 = 0.54. 

The calculations above were made for an ideal super- 
lattice, with no interfacial roughness. STE M 34 i 35 and 
HXR measurements (e.g., in Ref. 5) have shown that 
interface mixing in complex oxide layers is well approx- 
imated by a Gaussian roughness. STEM measurements 
on the LNO-LCO SL (data not shown) are consistent 
with this assumption. We now consider a Gaussian in- 
termixing distribution at the interfaces. Using as up- 
per and lower limits the surface roughness a = 0.53 nm 
(Section II. A) and 1/2 ML, because of the mixing of 
7 ML and 8 ML superperiods, we obtain an interface 
roughness Oi = 4.25 ± 1 A. Including roughness ef- 
fects in a scattering model suppresses the scattering in- 
tensity at L = 2 more than at L = 1 and increases 
the value for — p^t) calculated from the measure- 
ments. Specifically, A# (with roughness included) is 
Ar = V c 3(2Wc) 2 A « 1.49A and (p 3 + -p^ t ) R = 0.8. 

HXR measurements showed that the superlattice is 
made of both 8 ML and 7 ML superperiods (Section 
II. A); a similar calculation as the above for a 7 ML su- 
perperiod gives (p 3 ^ — p ut)n = 1-1. The average of the 
two structures is 0.95. These calculations show a valence 
change of 50 — 100 %, pointing to a considerable modu- 
lation of the Ni valence between the inner and the outer 
LNO layers, much larger than that expected from chem- 
ical bonding of 1 — 10 %. The large valence change is 
possible in our case because of the high doping of the 
LNO layers with interstitial oxygen. Choosing a larger 
valence for Ni ions inside the LNO layers is justified by 
the oxygen affinity of LNO; the Ni d^ z 2_ r 2 and d x 2_ y 2 
orbitals are emptied preferentially away from the inter- 
faces. The case of a valence change of 1 is represented 
schematically in Fig. 8b. 

The layer scattering factor was obtained above by in- 
terpolating between nominal Ni 2+ and Ni 3+ form fac- 
tors statistically with weights {pi}- The same results 
obtain for a scattering model describing a fractional va- 
lence common to all Ni ions in a MO2 layer, in which 
a linear dependence of the Ni ion form factor on dop- 
ing for intermediate valences between Ni 2+ and Ni 3+ is 
used. Which model better represents the reality depends 
on whether the oxygen interstitials dope only the imme- 
diate Ni neighbor ions or a large number of ions in the 
MO2 layers. Our experiment cannot distinguish between 
these two cases. 

The large Ni valence modulation would suggest a simi- 
lar effect in the valence of the oxygen ions, in addition to 
their layer-dependent density shown in Fig. 8b. One pos- 
sible indication of this effect is the feature seen at 533.5 
eV (Fig. 5). It is a state in the LNO layers since it aligns 
with a peak in TEY spectrum (Fig. 4). In contrast to the 
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LNO MCP state, it has a large intensity at L = 2, similar 
to the Ni edge scattering. It has also been observed for 
very high doping in the La 2 - x Sr x Cu04 system, which has 
no interstitial oxygen^ This suggests that it is an apical 
oxygen state hybridized with Ni orbitals, responding to 
the large Ni valence modulation. 

IV. CONCLUSION 

We observed a non-uniform distribution of holes in a 
LNO-LCO superlattice, with doped LNO layers and es- 
sentially undoped LCO layers. As for bulk LCO and 
LNO, the holes in the LCO layers are oriented in the 
a — b plane, while the holes in the LNO layers have a 
substantial c-axis component. The difference in the hole 
affinities of LCO and LNO layers modulates the valence 
and the orbital occupation of the Ni ions in the LNO 
layers, with the Ni d x 2_ y 2^ and d 3z 2_ r 2^ orbital states 
preferentially occupied near the interfaces. Since the den- 
sity of oxygen interstitials can be continuously modified 
with standard techniques, the Ni orbital occupation near 
the interfaces in 214 LNO-LCO superlattices is tunable. 
In contrast to the ferromagnetism of this structure, this 
orbital occupation is more robust, being stable up to at 



least 280 K. 

It would be useful to explore combinations of 214 ox- 
ides (manganites, titanates, etc.) in which the orbital 
occupations near the interfaces can be changed by selec- 
tive oxygen doping so as to leave only one orbital unoc- 
cupied, similar to superconducting cuprates. For super- 
conducting samples it might be desirable to avoid ferro- 
magnetism. On the other hand, it would be interesting 
to investigate 214 systems with magnetic orders with un- 
compensated spins at interfaces, for instance A-type anti- 
ferromagnets, for which the magnetic interactions across 
interfaces are sufficiently strong to correlate the spins and 
create a magnetically active interface! 37 ' 38 
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